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1 Introduction 



It is well-known that boundary conditions imposed upon quantum fields may lead to 
divergent expectation values for local observables. A simple example is a massless, 
minimally coupled scalar field ip(t,x) which vanishes on the z = plane; 



0. (1) 



One finds that the renormalized expectation values of both (p 2 and of the energy 
density T tt diverge as z — > 0. Specifically, 

(V 2 ) = 1 



16 7T 2 Z 2 



and 



(Ttt) — 2 4 • (3) 

16 TT Z Z 4 

(Units in which % = c = 1 will be used throughout this paper. The metric tensor 
is taken to be rj^ = diag(l, —1, —1, —1).) The stress energy tensor of the massless 
minimally coupled scalar field is given by 

Tp, = d^<p d u <p - ^ rj^ d a (p d a <p> . (4) 

Similar divergences occur in the expectation values of the squared electric or magnetic 
fields, (E 2 ) or (B 2 ) near a perfectly conducting plane, although in this case the local 
energy density remains finite. When the conducting boundary is curved, then the 
energy density diverges on the boundary 0. 

Furthermore, there is a puzzling discrepancy between the Casimir energy for a 
minimal scalar field computed as the renormalized expectation value of the Hamilto- 
nian and as a spatial integral of (Ttt). Consider the case of two parallel plates with 
separation L on which the field vanishes. If we first form the Hamiltonian operator 
H — J T tt d 3 x, the result is the same for both the minimal and conformally coupled 
fields. This follows from the fact that the stress tensor for the conformal field, 

1 1 r -, 

@^ = d^ip d u ip - - r)^ d a (pd a <p - - \d^(<p d v <p) + d„((p d^ip) - 2r)^ d a (<p d a <p)\ , (5) 

differs from that for the minimal field, Eq. (|]), by a total derivative term which 
integrates to zero. If we find the renormalized expectation value of H, the energy 
per unit area is found to be — 7r 2 /(1440 L 3 ). However, if we attempt to compute this 
energy per unit Jo (T tt ) dz, the result is divergent. This discrepancy has led 

some authors to postulate the existence of a singular surface energy density, which 
would render the latter expression finite and equal to the former result. (Note that 
the surface energy densities which are of concern here are distinct from the surface- 
area-dependent terms in the regularized Casimir energy which can arise in particular 
regularization methods [|], |f.) 
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In curved or topologically nontrivial spacetimes, it is also possible for the renor- 
malized expectation value of the stress tensor to diverge on particular boundaries. 
An example is the Boulware vacuum state in Schwarzschild spacetime, for which the 
stress tensor diverges on the event horizon |6|, [7J. This divergence is usually inter- 
preted as indicating that this is not a physically realizable state. Other examples 
of divergent stress tensors include Misner space, where the divergence occurs on the 
Cauchy horizon ||. In this and similar examples, one is tempted to resolve the 
problem by regarding the spacetime itself to be unphysical. Indeed, this philosophy 
is the basis of Hawking's Chronology Protection Conjecture pH |, which argues that 
closed timelike curves are prohibited by the effects of divergent energy densities which 
would otherwise appear on the chronology horizon (the boundary between a region 
containing closed timelike curves and one without such curves). 

One may understand why the imposition of a boundary condition such as Eq. (|l|) 
on a quantum field can result in infinities. In the case of {p 2 ), renormalization means 
taking the difference of the expectation value in the presence of the boundary and in 
its absence. Normally, this removes the infinite part and leaves a finite remainder. 
However, on the boundary the formal expectation value of (p 2 ) is finite, so the sub- 
traction results in an infinite difference. We can also understand why quantities such 
as (if 2 ) and (T a ) become infinite on the boundary. The field <p and its time deriva- 
tive p are conjugate variables which satisfy an uncertainty relation. If <p is precisely 
specified, p is completely indeterminate, and (p 2 ) and thus (T tt ) = \{p 2 + {Vp) 2 ) 
are infinite. A state in which p is precisely determined at a point has infinite energy 
density at that point for essentially the same reason that a position eigenstate in 
single particle quantum mechanics has infinite energy. 

In the case of material boundaries, such metal plates, infinite values of (E 2 ) or 
other observables are presumably avoided because such boundaries are not perfect 
reflectors at all frequencies. A metal plate is a good reflector of electromagnetic 
waves at frequencies below the plasma frequency, but becomes relatively transparent 
at higher frequencies. Such a high frequency cutoff seems not to be available when 
the "boundary" is a feature of the spacetime structure. 

The purpose of the present paper is to explore an alternative mechanism for 
introducing a cutoff which removes singular behavior on boundaries. This is to allow 
the position of the boundary to undergo quantum fluctuations. One might expect 
that such fluctuations will smear out the contributions of the high frequency modes 
without the need to introduce an explicit high frequency cutoff. 



2 ((p 2 ) near a Single Plate 

Let us consider a plane boundary located at z = q. If we impose the boundary 
condition on a massless quantized scalar field p that it vanish on this boundary, the 
appropriate two-point function may be constructed as an image sum. The result is 

(p(x)p(x')) = G(x, x') = G (x, x') + Gr(x, x') , (6) 
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where 

Go( *'*' ) = - 47rW-Ax2) (7) 
is the empty space two-point function, with At — t — t', Ax 2 = |x — x'| 2 , and 

Gr{ - X ' X>) = 47r 2 [At 2 - Ax 2 - Ay 2 - (z + z' - 2q) 2 ] ' (8) 

The full two-point function, G(x,x'), vanishes whenever z = q or z' = q. The renor- 
malized expectation value of ip 2 is given by the coincidence limit of the renormalized 
two-point function, Gr(x,x'), 

^> = G " ( '-'> = - 16,»(U)' <9) 

and is singular at z — q. 

We now wish to allow the position variable q to fluctuate. This will occur if the 
mirror is treated as a quantum object with a wavefunction ip(q), and hence a position 
probability distribution of 

f(q) = Mq)\ 2 , (10) 

where 

/•oo 

f(q)dq = l. (11) 



The average over position of a function H(q) becomes 

POO 

(H)= H(q)f(q)dq. (12) 



Thus to find the mean value of (p 2 , we need to calculate (G R ). This is most easily done 
by expressing Gr in a Fourier representation, and then averaging the g-dependence 
of the mode functions: 

(G R (x, x')) = -777^ Re / — e*'-^" 3 ^ e **«(*+*') ( e -«*.«) , (13) 

2(27r) ,:i J 

where k t and x t denote the components of k and x, respectively, in directions parallel 
to the plate. 

To proceed further, we must specify the probability distribution, f(q). A conve- 
nient choice is a Gaussian peaked about q — 0, 

/(?) = ^e"^ 2 , (14) 

V 71 

which leads to 
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with 



<« 2 > = i 



(16) 



This probability distribution is the appropriate one to describe, for example, a plate 
in the ground state of a harmonic potential. Note that Eq. fll5|) is equivalent to the 
result which one obtains when taking the vacuum expectation value of the complex 
exponential of a free quantum field. (See, for example, Eq. (8) of Ref. |§.) This is to 
be expected, as a free quantum field in the vacuum state is equivalent to an infinite 
collection of harmonic oscillators in their ground states. 
If we use Eq. (|15j) in Eq. (|T3|), set x' = x, and use 



d 6 k = 2tt 



duuj . 



for integrands independent of the azimuthal angle, the result is 

1 



(G R (t,t')) 



2(2n) 2 (t-t>] 



where 



Re 



A 



dk e i(2zk z + (t-t')\k z \) e -2k 2 z A< 



if) 



(17) 



;is) 



(19) 



is the root-mean-squared displacement of the mirror. This integral may be performed 
in terms of the error function $ to be 



{G R (t,t')) 



iV2n 



\(27T) 2 (t-t')A 



+ e 



-(2 Z +t-i') 2 /(8A 2 ) $ . 



-(2z-t+t') 2 /(8A 2 ) $ i 



.t-H + 2z y 



t-t' - 2z" 
2V2A t 



2V2A 



(20) 



(Here and at many other places in this paper, the calculations were performed with 
the aid of the symbolic algebra program MACSYMA.) This quantity is finite in the 
limit that t' —>■ t : 



(G 



V2; 



+ 



1 



32Vtt 3 A 3 V V^A J ' 16tt 2 A 5 

Note that this expression is real, as may be seen from the fact that 

2i 



-$>(-ix) 



—= / e u du. 

'IX Jo 



(21) 



(22) 



The quantity (Gr) is ((p 2 ) in the presence of position fluctuations, and is finite 
for all z. For large z, we have 



(Gr) 



3A 2 



16 tt 2 z 2 lQn 2 z 4 
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+ 



z > A, 



(23) 



<4> 2 >A 2 



0.006 



0.004 



0.002 



-0.002 




z/A 



Figure 1: 

The mean value of p> 2 near a mirror undergoing Gaussian position fluctuations is 
shown. Here A is the characteristic width of the probability distribution. For z/A 
large, (p> 2 ) is approximately given by Eq. @, but it is finite as z/A — > 0. 



thus recovering the usual form far from the mirror. As z — ► 0, 

1 



(Gr) 



16 7T 2 A 2 



(24) 



and is hence finite. For intermediate values of z, (Gr) may be computed numerically, 
and is depicted in Fig. 1. 



3 The Energy Density near a Single Plate 

3.1 Gaussian Fluctuations 

We now wish to apply the procedure used in the previous section to find the energy 
density in the presence of a single plate, whose position is undergoing fluctuations 
with a Gaussian probability distribution. This energy density may be expressed as 



(p) = X - Jim (d t d' t + V x • V x ,) (G R (x, x')) . 

x' - -• X 



(25) 
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A repetition of the procedure used for {Gr) leads to 



(p) 



1 i im Re f — (u 2 - k 2 ) e ikt <^-^ e^-V e ik *^ (e~ 2ik ^) . (26) 



2(2vr)3 



If we now employ the relation for Gaussian fluctuations, Eq. (|15|), and perform the 
integrations as before, we find 



(p) 



2i 



7T , d 2 

— hm 



(2vr) 3 A V 2 t'-*t dudv {u + v 



e -n 2 /(8A 2 ) $ 



U 



'2V2A 



2V2A 



)}} 



(27) 



where u = t — t' — 2z, and v = t — t' + 2z. Explicit evaluation of the last expression 
leads to 



(p) 



192 7r 2 A 7 



/ 2^z(z 2 -3A 2 )e- z2 ^ 2A2 U$ 



V / 2A / 



2A(z 2 - 2A S 



(28) 



Far from the mirror, the energy density is that calculated without fluctuations: 

1 



(p) ~ " 

and near the mirror it is finite 

(P) 



16 7T 2 Z 4 



z > A, 



(29) 



1 



48 7T 2 A 4 



z -> 0. 



(30) 



The energy density, (p) is given as a function of z in Fig. 2. 

The remaining components of the expectation value of the stress tensor, {T^) 
may be readily obtained. This must be a Lorentz tensor formed from the metric rj^ 
and n fl n u , where n M = (0, 0, 0, 1) is the unit normal vector to the mirror. Hence, 



{Tfj,u) = Ft(z) T)^ + F 2 (z) 



(31) 



where F\ and F 2 are scalar functions of z. We see immediately that the transverse 
components are just minus the energy density: 



(T XX ) = (Tyy) = -(p). 

Furthermore, the conservation law, <9 M = implies that 

d 



dz 



[F x {z)-F 2 {z)\=U. 



(32) 



(33) 



We define the renormalized stress tensor so that (T)^) — > as z — » 00, which implies 
Fi(z) = F 2 (z), and hence 



(Tp V ) = (p) {Vfiv + n^nv) . 



(34) 
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Figure 2: 

The mean energy density, (p), near a mirror undergoing Gaussian position fluctua- 
tions is shown. The integral of this function over all positive z vanishes. 



As a consequence, the pressure normal to the mirror vanishes: 

(T zz ) = 0. (35) 

Finally, we note that one may calculate the integral of the right-hand-side of 
Eq. fl28|) explicitly and verify that 

/ (p)dz = 0. (36) 
Jo 

This confirms that the boundary fluctuations remove the apparent discrepancy be- 
tween the Casimir energies of the minimal and conformal scalar fields. Note that 
although (p) is negative both at large distances and at the mean position of the 
mirror, it is positive in a region near z = 0, as illustrated in Fig. 2. The positive en- 
ergy region can be regarded as the concrete realization of the positive surface energy 
density conjectured by Kennedy, et al JJ. 

3.2 A General Probability Distribution 

Now we wish to generalize our discussion to an arbitrary probability distribution. 
For later use, we will momentarily assume that the distribution function / is peaked 



S 



symmetrically around an arbitrary value of (q), and write / = f(s), where s = q—{q)- 
The average of a complex exponential function of q then becomes 



(e iaq ) = e ia{g) f(a) 



(37) 



where / denotes the Fourier transform of /: 



/(«) 



f{s)ds. 



(38) 



We may now use Eq. (|26| ) to express the averaged energy density for an arbitrary, 
symmetric probability distribution as 



(p) 



lim zr—z Re 



2tt 2 v-*t dudv 



u + v Jo 



dk y 



f(-2k z )+e- ik * v f{2k z ) , (39) 



where, as before, u = t — t' — 2z, and v = t — t' + 2z. We next use Eq. (pq) to 
re-express (p) in terms of /, and employ the relation 



to find 



where 



(p) 



dxe 1 



1 , d 2 
hm 



7r 2 f-+t dudv 



+ n 6(a) 



F(u) + F(v) 
u + v 



F(u) 



ds 



(40) 



(41) 



(42) 



'-oo 2s + u 

and the last integral is understood to be a principal value. 

Thus given the probability distribution f(s), we need only calculate F (the Hilbert 
transform of /), and then evaluate the derivatives and limit in Eq. ([0]). It is of interest 
to apply this formalism to the case of a compactly-supported distribution. A simple 
example is 

315 



256 s; 



(s - s ) (s + s ) , -s < s < s 



(43) 



and f(s) = for |s| > Sq. This function is chosen so that / and its first three 
derivatives are continuous at s = ±sq. Equation (IT) now leads to 



(P) 



512 7T 2 s 9 



{(2205; 



3150s^ 3 + 945s^) 



hi(z + s ) — In \z — so| 



-4410s z 4 + 4830s^ 2 - 672s^} 



(44) 



This function is plotted in Fig. 3. Again, we have that at large distances from the 
mirror 



(p) 



1 

16 7T 2 Z 4 



+ 



z 3> s • 



(45) 




Figure 3: 

The mean energy density, (p), is illustrated for the compact probability distribution, 
Eq. (f43|), where so is the width of the probability distribution. Again, the integral of 
(p) over positive z vanishes. 

Note that (p) has a cusp at z = sq. If we had chosen a distribution function for which 
any of the first three derivatives are discontinuous at this point, (p) would become 
singular there. Similarly, we could smooth out the cusp by matching more derivatives 
at this point. We may verify directly that the total energy again vanishes: 

/ (p) dz = . (46) 



4 The Squared Electric Field near a Single Plate 

As noted in the Introduction, the squared electric and magnetic fields diverge in the 
presence of a perfectly reflecting plate with no position fluctuations. Specifically, one 
has 



so the energy density vanishes: 



(p) = i((E 2 ) + (B 2 ))=0. (48) 
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We may now calculate the squared electric field in the presence of a fluctuating 
boundary. The two point function for the photon field is 

D^(x, x') = (0\A^{x) A u (x')\0) . (49) 

In the presence of the reflecting plate, this may be expressed as 

D^(x, x') = D^(x - x') + D%(x, x') (50) 

where Dq V {x — x') is the two point function in the absence of the boundary, and 
the renormalized two point function, D R u (x, x'), is the correction introduced by the 
presence of the boundary. In a particular choice of gauge, we have that 

D^ u (x-x') =r]^G (x-x'), (51) 

and 

Dg(x, x') = + 2n fl n v )G R (x, x') . (52) 

Here Go and Gr are the scalar two point functions given in Eqs. (|7p and (|8]), respec- 
tively, and = (0, 0, 0, 1) is the unit vector normal to the plate. The renormalized 
field strength two point function can now be obtained by taking the four dimensional 
curl in x and in x' of Eq. (0). The electric field part of this function is 

(Ei(x) E^x')) = (F M (x) F 0j (x')) = dMMx) A 3 {x')) + d^A^x) A (x')) . (53) 

In the presence of boundary fluctuations, the mean squared electric field is now given 
by 

(E 2 ) = Jim (d t & t - V x • V x .) (G R (x, x')) = i<r;> , (54) 

x' — - X 

where is the trace of the minimal scalar field stress tensor. From Eqs. (^) and 
( |3~5|) , we find that 

(E 2 ) = -3(p). (55) 

Thus our explicit results for the scalar field energy density, (p), Eqs. (|28|) and (]44]) , 
also give us (E 2 ) for the Gaussian and compact probability distributions, respectively. 

It is of interest to note that the Casimir-Polder potential between a polarizable 
particle with a frequency-independent polarizability a and a conducting boundary is 
given by 

V(z) = ~a{E 2 ). (56) 

Thus Figs. 2 and 3 are also plots of 2V(z) / (3a) for the Gaussian and compact 
probability distributions, respectively. In both cases, there is a minimum in V(z) at 
a finite distance from the mean position of the boundary, at which the particle could 
apparently become trapped. This should probably not be taken too seriously. In the 
case of Gaussian fluctuations there is a nonzero probability to find the mirror to the 
right of the minimum of V(z). In the case of the compact probability distribution, 
Eq. (|43|) , the mirror may be found at any location to the left of the minimum of V(z). 
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5 The Energy Density between Two Plates 

Here we will address the problem of finding the mean energy density for a minimally 
coupled, massless scalar field between a pair of parallel plates. As before, the field is 
assumed to vanish on the plates, but their positions are allowed to fluctuate. First 
suppose that the plates are fixed at z = 77 and z — a, respectively. The two point 
function may be constructed as an image sum: 



G(x, x') 



1 



4tt 2 

00 

E 

m=— 00 



E 



1 



An 



[z — z' + 2m(cr — 77)]" 



At 2 -[z + z' + 2ma - 2(m + l)rf\< 



(57) 



where At 2 = At 2 — Ax^ = At 2 — Ax 2 — Ay 2 . The renormalized two point function 
is obtained by subtracting Gq(x,x'), which amounts to omitting the m = term in 
the first summation. The result may be written as 



G R (x,x') 



X 



2(2vr) : 



Re 



d 3 k 



UJ 



E < 

m=— 00 



ik z [z— z'+2m(<7— rj)] 



E < 

m=—oo 



ik z [z+z' +2ma— 2(m+l)rj] 



(5* 



where the prime on a summation denotes that the m = term is omitted. 

We will assume that the positions of both plates are described by the same prob- 
ability distribution, /. We then have that 



(e 



ik z [z—z'+2m(cr—r])] 



) 



drjf(rj) j daf(a) e u°.l*-*+*n(r-n)] 

ik z [z-z'+2m((a)-(r,))] p(2mk z ) , 



(59) 



and 



^«k z [z+z'+2mff-2(m+l) t) ]| __ e ik z [z+z' +2m(a)-2(m+l){ri)} f(2 m k z ) f (2(m + l)k. 



(60) 



The mean energy density is given by Eq. 
the above results, we find 



If we combine this expression with 



(p) = - 

idz 2 



lim Re 



in 



2 At^O 



dt 2 ^AtH- 



dk z e ilkzlAt e 2i|fcz|ma f 2 (2mk z ) + 

m=—oo 



vAt) J dkzeilkzlAt E e 2ilk > lima+z) f(2mk z ) f(2(m+ l)jfe. 



m=—oo 



,(61) 



where we have set (a) = and a = (77), so a is the mean separation between the 
plates. 
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We now employ Eq. PO] ) and the identity |LT 

1 



E 



m — a 



—Ti cot na . 



After the derivatives have been evaluated, our final result may be expressed as 

(P) = (P>1 + (P>2, 



where 



and 



(Ph 



71' 



1440 



(P>2 = " 



71 



48 J- 



ds 



ds 



drf(s)f(r) 



dr ■ 



mm 

(s + r + a 



,4 : 



2 sin 2 


n(z+r) 
_ s+r+a _ 


-3 


(s + r + a) 4 sin 4 


7r(z+r) 
s+r+a _ 



(62) 

(63) 
(64) 

(65) 



Let us first discuss the limit in which the position of both plates is precisely 
defined. In this case, we take f(s) = S(s) and obtain 



(Ph 



TV 



1440 a 4 



and 



2 2 sin 2 (&) - 3 



4 • 4 

a 4 sin 



(?) ' 



(66) 



(67) 



This is just the usual result 0; (p)i is now the energy density for a conformal scalar 
field, and (p)2 diverges on the boundaries. 

Now suppose that we take f(s) to be a function with a finite width, and whose 
first three derivatives are finite. The integrals in Eqs. (|64f) and ( |65"D contain poles 
in the ranges of integration, but the integrals are well-defined as principal value or 
generalized principal value integrals ||12||. That is, we use identities of the form 



ds 



m 



(s — a) 



ds 



Is — a) 



1 

6J-c 



ds 



f"(s) 



Is — a) 



ds 

oo S 



f"(s) 



(6£ 



Thus we see that both and (p)2 will be finite everywhere. 

We may now integrate the finite energy density on z to obtain the mean energy 
per unit area: 

E = 



where 
and 



(p) dz = E 1 + E 2 
Ei = a(p)i , 



/oo roc 
ds / drf(s)f(r)H(s,r) 
-oo J —oo 



(69) 

(70) 
(71) 
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and H(s,r) is denned by 



H(s,r) 



7i 

48 



cos 


7rr 


sin 3 


7r(z+r) 


— sin 3 


7rr 


COS 


7r(z+r) 


s+r+a 


s+r+a 


s+r+a 


s+r+a 


(s + r + a) 3 sin 3 


7rr 


sin 3 


ir(z+r) 




s+r+a 


s+r+a 



(72) 



In order to discuss the case of plates which are highly localized in position, we need 
to Taylor expand H(s, r) around s = r = 0: 



H(a,r) 



48tt 2 



7T 



\s + r) 



720a 4 



+ 



(73) 



Although the leading term in this expansion is singular at s = or r = 0, its 
contribution to E 2 vanishes because of the symmetry of the probability distribution, 
f(s) = /(— s). All subsequent terms in the expansion of H(s,r) vanish at s = r = 0. 
Thus we find that if we first form the total energy per unit area of the plates, and 
then take the limit in which f(s) — > 8(s), the result is the same as for the conformal 
scalar field: 



E 1 



71 



(74) 



1440 a 3 ' 

In both cases, we now find the same, negative Casimir energy. (Note that, in general, 
the sign of a Casimir energy is very difficult to predict in advance of an explicit 
calculation, and can depend upon both boundary conditions and the dimensionality 
of spacetime fll3~| .) 



6 Summary and Conclusions 

In the previous sections, we have seen that position fluctuations of a reflecting bound- 
ary are capable of removing divergences in the renormalized values of local observ- 
ables, such as (y? 2 ) and (T^). In the case of the massless, minimally coupled scalar 
field, such fluctuations also remove the discrepancy between the spatial integral of 
the renormalized energy density, / (T tt )d 3 x, and the renormalized expectation value 
of the Hamiltonian, (H). Position fluctuations are necessary if one is to treat the 
mirror as a quantum mechanical object. 

Of course, for real mirrors the mass is likely to be so large that the position 
uncertainty A is very small. In this case, the cutoff in reflectivity due to dispersion 
will normally be the dominant effect. Dielectric materials become transparent to 
electromagnetic radiation at wavelengths shorter than about the plasma wavelength, 
A p . So long as the position uncertainty is small compared to this length, A <^ 
X p , dispersive effects are dominant, and the position fluctuations may be ignored. 
However, if one could arrange to prepare a mirror in a quantum state in which A > \ p , 
the position fluctuation effects discussed in this paper would become dominant. 

In the area of gravitational physics, the situation is rather different. Here it is 
also possible to have horizons which act as boundaries for quantized fields, and for 
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the renormalized expectation value of the stress tensor to diverge on the horizon. As 
discussed in the Introduction, examples include the event horizon of Schwarzschild 
spacetime in the Boulware vacuum, the Cauchy horizon in Misner space, and possibly 
the chronology horizon in a spacetime containing closed timelike curves. Now there 
is no natural cutoff at high frequencies, and in fact higher frequency modes tend 
to couple more strongly to gravity by virtue of their larger energy. It is of course 
possible that a more complete quantum theory of gravity will introduce an effective 
cutoff at the Planck scale. At the present, any discussion of Planck scale physics must 
be regarded as highly speculative. Nonetheless, position fluctuations of the horizon 
would seem to provide a possible way to avoid divergent stress tensors. It is plausible 
that the location of a spacetime horizon undergoes position fluctuations due either to 
the quantum nature of gravity ("active fluctuations") [[14]], or to fluctuations of the 



stress tensor of quantum matter fields ("passive fluctuations") [jn|. This is a topic 
requiring further study. 
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